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§J Reversible dissipation in Euler and Navier Stokes equations. 

In [Ga6] one finds an analysis leading to a conjecture on the equivalence be- 
tween the irreversible NS equation and a reversible equation that was called 
GNS ("gaussian NS equation"). 

The ideas discussed in [Ga6] are, however, much more general and it is worth 
pointing out some other applications, as well as several possible tests that seem 
under reach of present day experimental (numerical and "real") techniques. 

I shall focus on fluid mechanics problems considering a fluid that: 

(1) is enclosed in a periodic box with side L, possibly with a few disks 
(" obstacles" ) removed so that no infinite straight path can be found in SI that 
avoids the obstacles, 

(2) is incompressible with density p. 

I shall consider four distinct evolution equations for this fluid, all of dissipative 
nature. 

u+u -d u = — —dp + g + i/Am , d ■ u = NS 
P 

u+U 9 u = —dp+ g +/3Au, d-u=0 GNS 

\ ~ " ~ ~ (1.D 
u+u ■ d u — — dp + g — xik-, d ■ u = ED 
P 

u+u ■ d u — dp + g — au, d ■ u — GED 

P 

In the case f2 contains obstacles a "no friction" boundary condition will be 
imposed on dfl, i.e. u ■ n — if n is the normal to dfl. The first equation is 
the well known Navier Stokes equation with v being the viscosity. 
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The second equation, introduced in [Ga6] and called GNS, has a multiplier (3 
defined so that the total vorticity nL 3 = p J lu 2 dx, with u = d A u being the 
vorticity, is a constant of motion; this means that: 

/ A g ■ lj + uj_ ■ (lu ■ d u )) dx 
= J^dx " ™ 

The third equation will be called the Euler dissipative equation, ED: it repre- 
sents a non viscous ideal fluid moving in a " sticky background" : x is a " sticky" 
viscosity. The model is, as far as I know, not a good model for any physical 
situation, but it is interesting to consider it mainly for comparison purposes. 

The fourth equation will be called GED equation, gaussian dissipative Euler 
equation and here a is a multiplier defined so that the total (kinetic) energy 
eh 3 = | J u 2 dx is a constant of motion in spite of the action of the force g ; 
this means that a is given by: 

f q ■ u 

*(u) = L =^- (1-3) 

A similar equation, with a constraint on the energy contained in each "mo- 
mentum shell" to be constant was considered in [ZJ], which is the first paper 
in which the idea of a reversible Navicr Stokes equation is advanced and stud- 
ied. The energy content of each "momentum shell" was fixed to be the value 
predicted by Kolmogorov theory, [LL] . 

Note that both the GED and the GNS equations have a symmetry in u , 
so that they are reversible in the sense that, if Vt is the flow describing the 
equation solution (so that t — > VtU = u(t) is the solution with initial data u), 
then the transformation i : u — > — u anti-commutes with the time evolution 
V t : 

iV t = V_ t i (1.4) 

We shall avoid (as it is, unfortunately, always the case in the current litera- 
ture) considering the problem of proving the global existence of solution to the 
equations (1.1) (the problem is in fact open, see [Gal]) and we shall consider 
the truncated equations with momentum cut off K. 

The truncation will be performed on a suitable orthonormal basis for the 
divergenceless fields in £1: given the boundary conditions we consider it natural 
to use the basis generated by the minimax principle applied to the Dirichlct 
quadratic form J Q (d u) 2 dx defined on the space of the C°°(f2) divergenceless 

fields u with u ■ n = on dfl. The basis fields Uj will verify: Au^ = 
—EjUj + djp, for a suitable multiplier fij, with UpHj G C°° and Ej are 
eigenvalues) . 

For instance in the case of no obstacles let the u =X)fc^o 7 k e% ~'~ ^ e ^ e 
velocity field represented in Fourier series with 7 = ~ and k ■ 7 =0 

(incomprcssibility condition); here k has components that are integer multiples 
of the "lowest momentum" ko = ^j-- Then consider the equation: 
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ik=-#(k)l k - i (lfc t ' ^2)Hfe 7 i2 + 5^ (1-5) 

where the k 's take only the values < | fc | < if for some momentum cut-off 
K > and II ^ is the projection on the plane orthogonal to fc. This is an 
equation that defines a " truncation on the momentum sphere with radius K of 
the equations (1.1)" if: 

$(k) = -vk 2 NS case 
0U) = -/?fc 2 GNS case 

(1.6) 

0(k) = -X ED case 
■d(k) = -a GED case 

For simplicity we may suppose, in this no obstacles cases, that the mode fc = 
is absent, i.e. 7 = 0: this can be done if, as we suppose, the external force 

g does not have a zero mode component {i.e. if it has zero average). 
In order that the resulting cut-off equations be physically acceptable, and 
supposing that g =/= only for |fc| ~ fco, one shall have to fix K large. For 
instance in the NS case it should be much larger than the Kolmogorov scale 
K = (^z/ -2 ) 1 / 4 , where vr\ is the average dissipation of the solutions to (1.6) 
with K = +00 (determined on the basis of heuristic dimensional considerations 
by 77 - \g\ 2 L 2 v- 2 ): see [LL]. 

We shall use the same cut off for the other equations with fc replaced by the 
basis label j and | fc | replaced by \/~Ej, which is certainly a natural choice for 
the GNS equation. 

For the ED and GNE equations the choice of K should be made by developing 
a theory analogous to Kolmogorov's theory. We only attempt a preliminary 
analysis in §6 as the latter equations have no physical interpretation (of which 
I am aware) and they are used here only for the purpose of illustrating some 
interesting mechanisms and theories. Below we always refer to the truncated 
equations, unless otherwise stated. 

It is easy, in the no obstacles cases, to express the coefficients a, (3 for the cut 
off equations: 



2^0< k\<K 9 u " 7 k 

a=^^ P = P> + Pe (1.7) 

l^0<k\<K X k 

_ Efc^o k 2 g k ■ Tfc 
/e Ekk 2 h,\ 2 

a - J E il+i2+ fc 3= o klii ■ k 2 )(i k2 ■ i k ) 

Pi = ^— — — =2 - 

Ekk 2 \ lk \ 2 

where the fc 's take only the values < \ k\ < K for some momentum cut-off 
K > and life is the orthogonal projection on the plane perpendicular to fc. 
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The cases in which the region Q contains obstacles is very similar but we 
cannot write simple expressions for the basis fields and therefore the equations, 
although formally very similar to (1.5), (1.7) cannot be written very explicitly. 

The solutions of the equations (1.5), or to the corresponding ones in the obasta- 
cles cases, will be denoted V t u ' ns u , V^' 9ns u , V*' ed u , V^' ged u when the initial 
datum is u . Or in general: 

V^u, £= (v,ns), (rj,gns), (x,ed), (e,ged) (1.8) 

Keeping the forcing g constant we shall admit that for each equation, i. e. for 
each choice of £, there is a unique stationary distribution fj,^ describing the 
statistics of all initial data u that are randomly chosen with a " Liouvillc dis- 
tribution", i.e. (in the no obstacle cases, to fix the ideas) with a distribution 
Ho(d"f) proportional to the volume measure 5^ Yi \ k \ <k ^7 & ' w ^ ere the delta 
function is present only in the case of the reversible equations and fixes the 
constants of motion to the value prescribed by the first label in £. 

This means that given "any observable" F on the phase space T (of the 
velocity fields with momentum cut-off K) it is: 

lim i f F{vh )dt= I F( 7 ')^(d 7 ') d ={F) ( (1.9) 

for all choices of 7 except a set of zero Liouville measure. The distribution 
li ns will be called the SRB distribution for the Eq. (1.5), (1.6) (see the "zero-th 
law" in [UF],[GC2]). 

A particular role will be plaid by the averages (e)|, (77) $ as well as by the 
averages of (a){, (/3)| and of the entropy production rate cr(-f) that is defined 
by the divergence of the r.h.s. of the cut off equations. 

Consider explicitly only the no obstacles case: if Dk is the number of modes 
k with < I k I < K then the number of (independent) components of { 7 } 

is 2Dk and, see (1.5), setting 2Dk — E|/s|<ir ^ — 2 (which in the case with 
obstacles become 2Dk = J2^/W < k V^i))' 0110 nn d s that a is given by: 



a =2Dkv £ = {v, ns) 

o-=2D K (3-^ e -^ i ^=(t],gns) 

o=2D KX Z={x,ed) 

a =2Dkoi — a £ = (e, ged) 



(1.10) 



where f3 i ,f3 e are suitably defined, e.g. in the no obstacles cases: 

E,fc 2 liJ 2 (£ 4 * a l7 4 l a ) a 

so that a ~ 2DkP for £ — (V: 9 ns ) an d & — 2Dk& for £ = (e, ged). 
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The equivalence of dynamical ensembles conjecture, [Ga6], is the following: 1 

Conjecture NS: The statistics jJL v , ns , /J-q,gns °f the NS equations and of the 
GNS equations respectively are equivalent in the limit of large Reynolds number 
provided the parameters r\ and v are so related that (o~) v , ns = {o~) 7hgns . 

Here equivalent means that the ratios of the averages of the same observables 
with respect to the two distributions approaches 1 as R — > oo. The Reynolds 
number is defined here to be R — e 1 / 3 ^ 4 / 3 ^ -1 . 

A corresponding conjecture can be formulated for the ED and GED equations: 

Conjecture ED: The statistics H x , e d, f-e^ed of the ED equations and of the 
GED equations respectively are equivalent in the limit of large Reynolds number 
provided the parameters e and \ are so related that (cr) x , e d — (o-)c, ge d- 

The above stated conjectures are closely analogous to the familiar statements 
on the equivalence of thermodynamic ensembles, with the thermodynamic limit 
replaced by the limit R — > oo of infinite Reynolds number. They can be 
substantially weakened for the purposes of possible applications. 

The idea of non equilibrium ensembles and their possible equivalence is not 
really new: the recent literature contains many hints in such direction. The 
clearest is perhaps [ZJ]. See also [Ga2] (§4) and [MR]. 

On heuristic grounds, the conjectures would be justified if one did accept 
that the entropy creation rate reaches its average on a time scale that is fast 
compared to the hydrodynamical scales. The coefficients a ~ (2Dk)^ 1 <J, and 
ft ~ (2Dk sec (1.10), would be confused with their time averages {a) e ^ gne 
or (/3) v , gns and identified with the viscosity constants v or \- 

In this way the GNS and the NS equations would be equivalently good: both 
being the macroscopic manifestation of two equivalent microscopic dissipation 
mechanisms: one explicitly specified by the Gaussian constraint of constant 
dissipation and the other with dissipation unspecified a priori but phenomeno- 
logically modeled by a constant viscosity. Likewise one can view the GED and 
the ED equations as macroscopically equivalent: one with constant energy and 
the other with constant sticky viscosity \- 

The interest of the above conjectures is that the same physical system in 
which irreversible dissipation occurs (the NS or ED equations) can be described 
equivalently by a reversibly dissipative system (the GNS or GED equations). 

For reversible systems a general principle, the chaotic hypothesis, can be rea- 
sonably assumed to hold and to imply consequences that seem to be non trivial, 
see [GC1], [GC2], [Ga3], [Ga4],[Ga5]. 

The next section is devoted to a quick discussion of some of the established 
consequences of the principle and §3%§6 will deal with heuristic ideas and 
with describing a possible scenario for the phenomenology of the equations 
(1.1). The scenario will be developed without any pretension of rigor and it 



... e tanto nuova e, nclla prima apprensione, remota dal verisimile, che quando non si 
avesse modo di dilucidarla c rcndcrla piii chiara che'l Sole meglio sarebbe il tacerla che'l 
pronunziarla; pero, gia che me la son lasciata scappare di bocca..., [G], p. 231. 
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will present what will appear as the simplest among many other possibilities. It 
leads (implicitly) to several possible experimental tests of the chaotic hypothesis 
and of the other ideas involved in its development: the tests can also be viewed, 
independently just as interesting experiments proposals. 
In the following we shall always consider the NS equations and the GNS equa- 
tions with parameters fixed so that \i v , ns and fi v ,gns are equivalent by the con- 
jecture NS, and likewise we shall always consider the ED and GED equations 
with parameters fixed so that H x ,ed and fi £ ,ged are equivalent by the conjecture 
ED. 

$2 The fluctuation theorems. 

In reference [GC1],[GC2] the chaotic hypothesis was presented as a reformu- 
lation of an older principle due to Ruche, [Rl]. It gave us the possibility of 
some quantitative parameterless "predictions", in various cases, see also [Ga3], 
[Ga4], [Ga5]. The hypothesis is: 

Chaotic hypothesis: A chaotic many particle system or fluid in a stationary 
state can be regarded, for the purpose of computing macroscopic properties, as 
a smooth dynamical system with a transitive Axiom A global attractor. In re- 
versible systems it can be regarded, for the same purposes, as a smooth transitive 
Anosov system. 

The main result of [GC1] is the fluctuation theorem that gives a property of the 
variable p = p( 7 ) defined in terms of the contraction rate (To of the attractor 
surface elements by: 

1 r' 2 

- / a (V tl )dt = (a )+p (2.1) 

T J-r/2 

which can be regarded as a random variable with the distribution -K T {p)dp that 
it inherits from /i,^ ns or ^ e ,ged] here (a ) + is the average of a with respect to 
the distribution n v ,gns or n e . ged . 

Note that a should not be confused with a, (1.8): thinking of the attractor 
as a smooth surface cr is the contraction rate of its surface elements, which is 
different from the contraction rate a of the phase space volume elements, sec 
§4%§6. 

If the conjectures of §1 are accepted (cr )+ is also the /i„, ns or [i x . fi d average 
of do or at least tends to it as R — > 00 . 

If (<7o)+ > 0, see [R2] for a discussion of the conditions for this inequality 
("Ruelle's H-theorem" ) , and if ((p) = lim^oo ^ log7r r (p) then the fluctuation 
theorem of [GC1] gives the following large deviation relation, see also [Ga3], 
[Ge], for the equations GNS and GED: 

COO - CHO = 1; forallp (2 2) 

{0-0} +p 

which, in the case of noncquilibrium statistical mechanics, has been interpreted 
as an extension of the fluctuation dissipation theorem to large forcing fields, 
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[Ga5]. Here "for all" p means for all possible values of p (which is in general a 
bounded quantity). 

The fluctuation theorem (2.2) says that the distribution of p is multifractal, 
not surprisingly since £(p) can be regarded as a kind of generalized sum of 
Lyapunov exponents in the sense of [FP], [BJPV], and the odd part of ((p) is 
linear. 

A more general fluctuation theorem concerns the joint distribution of the 
variable p and of any other variable 5 = 9(7) that is similarly defined in terms 
of an observable Q which is odd under the time reversal operation that is defined 
on the attractor, i.e. : 

- f ' dtQ(vh) = (Q)+q (2-3) 

T J-t/2 

If n T (p,q) denotes the joint probability density of the observables p, q and if 
C(P) 9) = um T— >oo log tt t (p, q) then it follows from the chaotic hypothesis that 
the distributions for p, q, with respect to the statistics /j,^, £ = (f3,gns) or 
£ = (a,ged), verify: 

C(p, q) - C(-p, -q) , c „ /„ 

— = 1, for all p, q 2.4 

which, in the case of nonequilibrium statistical mechanics, has been interpreted 
as an extension of the Onsager's reciprocity to large forcing fields, [Ga5]. 
In the case of the equations (1.1), second and fourth, the above relations are 
applicable when the motions are chaotic, i.e. have at least one positive Lya- 
punov exponent (which should happen as soon as R is large enough, excepted 
possibly very special cases, see [Ma]); and it gives an interesting parameterless 
prediction if the contraction rate cto can be related to the contraction rate a of 
the GNS equations. 

If the conjecture in §1 holds in the strong sense of complete asymptotic equiv- 
alence between the "ensembles" fi-q.gns and (i^ns (or fJ- e .ged and /i x , e d) then 
(2.2) should also hold for other models of the viscous stationary states, like the 
one given by the classical NS equation in particular, note the quantities a in 
(1.10) are still fluctuating variables even when the evolution considered is given 
by the NS or ED equation. A check of this property might be under reach of 
present day experimental techniques. 

A check of (2.2) or (2.4) might be more accessible in the case of fluid systems 
because they often show chaotic motions with relatively few degrees of freedom 
so that the large fluctuations, that must occur in order to make possible testing 
the fluctuation theorems, are more likely occur and be observable. 

But this also means that the attractors have a very small size, compared to 
that of phase space: hence a problem in the interpretation of (2.2) is the fact 
that the quantity <To that appears in the theorem does coincide with the easily 
determined (see (1.10)) contraction of volume in the phase space only if the 
attractor is dense in phase space. 
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This is a property that can be expected in the case of noncquilibrium statistical 
mechanics under small or moderately strong forcing, see [BGG] for a discussion 
of this point, but it cannot be expected at large forcing or in fluid systems. In 
the latter case one should use the contraction rate of the volume elements on 
the attractor, see [BGG], [BG]. In the GNS systems it is likely that the property 
never really holds. 

This might render (2.2) quite useless as it is usually unrealistic to hope to 
determine the attractor equation with accuracy sufficient to compute its area 
contraction per unit time (assuming, as the chaotic hypothesis implies, that 
the attractor can be regarded as a smooth surface). 

Nevertheless in nonequilibrium statistical mechanics further analysis is possi- 
ble, based on an important symmetry of the Lyapunov exponents spectrum, 
and one can relate easily the area contraction on the surface defined by the 
attractor and the total phase space volume contraction. Hence one can try to 
push further the analysis of [BGG], [BG] to sec if one can say something also 
in the case of the GNS equations. 

§5 A scenario for experimental checks of the fluctuation theorem. The Lya- 
punov spectrum of ED and GED equations. Pairing rule and Axiom C. 

What follows is a very heuristic analysis aimed at giving an argument for the 
explicit form that the fluctuation theorem (2.2) will take in the case of the 
GED equations: and therefore by the equivalence conjecture also for the ED 
equations. No pretention of mathematical rigor is present and the idea is to 
illustrate the simplest possible scenario that I consider possible. The interest 
is (apart from the subjective feeling of a certain beauty) that the discussion 
suggests experiments and checks that have intrinsic interest and that do not 
seem to have yet been considered in the literature. 

We consider first the case of (1.1) in a domain Q, with obstacles: in spite of the 
appearances this is an easier case because in this case we can imagine forcing the 
system with a locally conservative force which is not globally conservative, like a 
field roughly parallel to one axis and tangent to the obstacles (one can imagine 
a uniformly charged fluid under an electromotive constant electric field). 

Note that in order to have a non trivial forcing the forcing field must be non 
globally conservative: otherwise its effect would be just that of altering the 
pressure. 

The Euler equations can, in general, be regarded as hamiltonian equations 
for a system whose configurations are the diffeomorphisms of the box Q (in 
our case a torus with, possibly, a few holes) containing the fluid: they are not 
directly in hamiltonian form in the same sense as the (closely analogous) Euler 
equations for a rigid body with a fixed point are not immediately hamiltonian 
(e.g. they involve half the number of actual equations of motion). 

In this way the GNS or GED equations can be regarded as hamiltonian equa- 
tions (approximately so, because the cut-off K destroys this property) modified 
by the action of a non conservative force g and by the gaussian constraint that 
the total vorticity or the total energy are constants. 
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Of course we exploit the "slip" (i.e. no friction) boundary conditions in order 
to be able to conclude the hamiltonian nature of the Eulcr equations. 

The phase space will then consist of a space larger than the above T, see (1.9): 
its points (u, 6) will be (cut-off) velocity fields and (cut-off) displacement 
fields describing the positions of the fluid particles with respect to a reference 
configuration. We call this the "full phase space" of the equations (1.1). 

The equations for the displacements will be in all models (1.1): 

S = u(6,t), S(x,0)=S (x) (3.1) 

which, once u is known from (1.3), (1.2), permit us to compute the physical 
fluid flow and the positions S(x,t) of the fluid particles that at time were 
at the points 5 (x), away from the reference configuration position x G f2. 
In the case in which the u verify truncated equations also (3.1) have to be 
truncated, for instance by replacing each e 1 -'-^-^ in u(5_,t) by its truncated 
Fourier expansion. 

The system motions (describing velocity and displacement fields) can be re- 
garded as motions with 2Dk degrees of freedom where, for instance in the no 
obstacles case, Dk is the number of non zero modes k with < \k \ < K (be- 
cause each 7 has two complex components but 1 _ k =~ k )- This means 
that ADk coordinates are necessary to describe the motion. 

Hence there are ADk Lyapunov exponents, 2Dk from the velocity equations 
(1.2) and 2D K from the displacements equations (3.1). 

In view of the equivalence conjectures we study the equations GNS and GED 
when convenient and the NS or ED when convenient. 

Out of the ADk exponents one has to extract, in the GNS or GED cases, 
one exponent that is trivially because of the conservation of the dissipation 
rate and one exponent that is trivially zero and corresponds to the vector field 
given by the r.h.s. of the GNS equation. Furthermore in the GNS or GED cases 
two more vanishing Lyapunov exponents are associated with other constants 
of motion. 1 

The other 2N = AD K - A exponents, or in the ED, NS cases all the 2N = 
ADk exponents, can be ordered in two groups the first containing the first 
N exponents in decreasing order and the second the remaining N ones in 
increasing order. 

The exponents of the first group are denoted A^~, j = 1, . . . , N and the ones 
in the second group are denoted Xj , j = 1, . . . , N. We call the two exponents 
(A+, Xj) a pair. 

We consider first in detail the ED and GED equations. In the above context 
it seems reasonable that in the full phase space of the GED and ED equations 
a pairing rule holds: 

K + K 

3 1 = const (3.2) 



1 In the case of the GNS equations the helicity J u ■ u dx is a constant of the motion and 
such is J 8 2 dx for the displacement equations. 
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at least when the forcing is locally conservative as we suppose from now on 
unless otherwise stated. The value of the constant will be called the "pairing 
leveF or "pairing constant', which must be ^(a) + , see (1.10). 

The pairing rule, in fact, formally holds in the present ED case. One can easily 
adapt the proof in [Dr]: this is discussed in the Appendix Al. 

The rule then holds also for GED as a consequence of the conjecture ED. A 
direct proof can probably be made along the lines of the work [DM] . In fact the 
constraint imposed by the definition of the multiplier a, (1.7), is a constraint 
of the type called isokinetic in [DM] and their proof seems to apply "without 
change" , although I did not check the details (the appendix Al should give the 
background for such an analysis). 

In the cases in which (3.2) has been proved, [DM], [Dr], it holds also in a far 
stronger sense: the local Lyapunov exponents, 2 of which the Lyapunov expo- 
nents are the averages, are paired as in (3.2) to a constant that is j independent 
but, of course, is dependent on the point in phase space. We call this the strong 
pairing rule. See the final comments. 

Note that the Lyapunov exponents of the full system can also be easily divided 
into velocity exponents, i.e. the ones of the GED or ED equations, and the dis- 
placement exponents, i.e. the others (which cannot be measured from the GED 
or ED evolution alone but require also (3.1)). In fact if we denote symbolically 
by (x, y) the pair ( u , 6 ) then the jacobian matrix of the equations is described 



For a further classification of the exponents we shall think that the Lyapunov 
exponents are divided into three classes that we call viscous, inertial and slow. 
The following scenario will be again summarized and enriched in the figure in 



(1) The slow exponents ("slow pairs") consist of M pairs of exponents the 
largest of which is < and it is a velocity exponent corresponding to slow 
motions of the velocity field, while the other (necessarily < 0) exponent of 
the pair is a displacement exponent and corresponds to a fast approach to the 
stationary state of some of the displacement variables. 

(2) The viscous exponents ("viscous pairs") consist of V negative velocity ex- 
ponents describing the fast approach to the stationary state of the viscous 
degrees of freedom of the velocity field: their paired positive exponents are 
displacement exponents associated with chaotic motions of the displacement 
variables. 

(3) The remaining 2P = N — M — V pairs (" inertial pairs" ) have one > and 
one < Lyapunov exponents: P of the pairs are pairs of velocity exponents 
and the P other pairs are displacement exponents. The P pairs of velocity 
exponents are the only pairs of exponents of the equations for the velocity 
field that contain one positive and one negative element: they describe the 



2 i.e. the logarithms of the eigenvalues of \J jf Jt , if Jt is the local jacobian matrix of the 
evolution operator Vt , other than those relative to the directions of the flow or of the imposed 
constraints or of other constants of motion 




where A, B, C are operators. 



§4. 
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gross characteristics of the chaotic motion on the attractor. It follows that the 
three types of exponents can in principle be uniquely identified among the N 
exponents of the velocity field equations, see also below. 

The existence of a certain number denoted P above of pairs of exponents, for 
the velocity field evolution, that are pairs of exponents of opposite sign does 
not follow simply from the fact that we are collecting togheter pairs containing 
a > exponent. In principle the > exponents of the velocity field could be 
paired with negative displacement exponents. We think that it is natural that 
the > Lyapunov exponents for the velocity field are paired with < exponents 
of the velocity field because we associate such pairs with the motions on the 
attractor. Since the GED equations are reversible it follows from [BG] that 
if the motions are also supposed to verify a geometric property called in [BG] 
Axiom C property (a simple extension of the paradigm of turbulent behavior, 
see [Rl] , that is the Axiom A property) then there must be an equal number P 
of positive and negative exponents for the restriction of the GED equations to 
their attractor. It seems therefore natural to think that they form P pairs. 
The equality of the number of > and < exponents for the motion on the 
attractor for the velocity fields is due to the existence, in reversible Axiom 
C systems, of a local time reversal map i* that transforms the attractor into 
itself anti-commuting with the time evolution, even when (and this is the rule 
in fluid dynamics) the attractor itself is not time reversal invariant: see [BG]. 
We proceed under the assumption that the Axiom C property is verified: for a 
complete discussion of the property we must unfortunately refer to [BG]. 
In Axiom C systems the time reversal symmetry "cannot be lost": when it is 
spontaneously broken (because the attractor is not time reversal invariant) it is 
replaced by a "weaker" symmetry, good enough to make "effectively reversible" 
all the motions on the attractor, a relation similar to the one in fundamental 
Physics between T and TCP (the latter being the "real" time reversal as the 
first is not a symmetry of the world we see). 

(4) The other P pairs should consist of displacement exponents exhibiting a 
rather symmetric behavior with respect to that of the GED exponents. Below 
we are going to suggest that very similar properties hold for the NS and GNS 
equations: in that case this further appealing symmetry seems compatible with 
(and in fact it was suggested by) the data on the velocity Lyapunov spectrum 
for models (" GOY shell models" ) whose behavior is " believed" to be related to 
NS equations: see [BJPV], figure in p. 71, taking into account that the pairing 
level in such data is very small because the viscosity is very small. 
In the above scenario the existence of the other P pairs of displacement expo- 
nents is assumed in order to make the total count of the number of exponents 
correct and is not based on evidence of any other kind. The displacement ex- 
ponents do not seem to have been considered in the literature, at least not in 
conjunction with the velocity field exponents (not surprisingly in view of the 
difficulty of the measurements) . 
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Thus if we measure the Lyapunov exponents for the GED equations alone we 
expect to find P pairs of opposite sign exponents paired at the value -^(a}+ 
for some P < N. 

It seems reasonable that the P pairs of displacement exponents coincide with 
the P pairs of inertial exponents for the velocity field equations: but this is not 
really necessary in order that an unambiguous identification of the three type 
of exponents be possible. They are already identified by the above properties. 

However if the P pairs of velocity inertial exponents and the P pairs of dis- 
placement inertial exponents do coincide we see that, by the pairing rule, the 
knowledge of the Lyapunov spectrum for the velocity equations implies that all 
the displacement exponents are knwon as well: no need to compute them. 

§^ Predictions of the fluctuation theorem for ED and GED. 

With the scenario developed in §3 we reconsider the fluctuation theorem and 
note that it is easy to check, by evaluating the divergence of the r.h.s. of the 
equations (1.2), (3.1), that the volume in the full (27V dimensional) phase space 
contracts at the same rate a at which the volume in velocity space does. Fluid 
incompressibility, and absence of displacement variables in the equations for 
the velocity field, imply this property. 

Furthermore if the strong pairing rule is assumed the total volume contraction 
in the full phase space, including the displacement variables, will be (7(7) and 
it will be related to the contraction o~o ( 7 ) of the area on the attractor surface by 
00(7) = ■jft°~(l), see [BGG] where the same mechanism was first exploited. 
This gives proportionality between the " apparent" contraction rate a and the 
"true" contraction rate 00 on the attractor for the GNS equations. 

As discussed at the end of §2 the fluctuation theorem holds for the fluctuations 
ofao so that the fluctuations of a will verify (2.2) but with a r.h.s in which 1 is 
replaced by where P is the number of pairs of Lyapunov exponents for the 
GED equations with one positive element. 

If the number of degrees of freedom is increased by increasing K one should 
expect, therefore, that the constant 757(0")+ approaches P(a}+ because the 
number of "true exponents" (i.e. inertial exponents) should not change as soon 
as K is so large that the motion is well described by the truncated equations: 
in fact it is believed that the attractor dimension does not depend on the 
truncation scale K (as long as it is large enough). 

Since the conjecture of §2 implies (a)+ = x, the constant should approach P\, 
at least if R (i.e. the Reynolds number) is large. The fluctuation theorem will 
thus take the form: 



«P>-«-P) =px (41) 

if the variable p is defined as in (2.1) but with the measurable a replacing the 
a priori difficult to measure a$, and £(p) is the limit of t -1 log7r T (p) with ir T (p) 
being the ti e ,ged distribution of p. The number P is accessible by measurements 
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performed only on the GED equations and not involving the displacement vari- 
ables (being the number of positive Lyapunov exponents of the GED equations) . 
The above analysis is somewhat conjectural but experiments, at least numerical 
ones, are possible to check the picture; e.g. one could attempt at: 

(1) checking the just derived slope P\ or 

(2) checking the following picture, representing the above classification of the 
exponents: 



pairing 

level 



A 



X M 



A 



X j-, 



A 



The continuous line in the first graph gives the value (< 0) of the j-th (among 
M) slow Lyapunov exponent (as a function of jj) of the GED equations; the 
dashed line is the graph of the paired exponents (of the displacement equations) 
and the intermediate line is the pairing constant. The exponents are defined 
only for x — jj but the graphs give, instead, continuous (or dashed) lines for 
visual aid. 

The second graph gives the values of the j-th (among P) pair of inertial 
exponents of the GED equation (one positive and one negative per pair): here 
too we use the continuous curves even though the number of such exponents 
will ususally be much smaller than the total number and therefore a discrete 
representation would be more appropriate. 

The negative curve in the third graph is the graph of the j-th viscous exponent 
(out of V) of the GED equation; the corresponding positive curve (dashed) 
is the curve of the companion exponents which correspond to displacement 
exponents. A fourth graph giving the other P displacement exponents (in 
pairs of one negative and one positive) would be qualitatively equal to the 
second graph (with the curves dashed for consistence of notation) . 

The graphs are not experimental data: they are just sketches illustrating the 
"simplest" picture that I considered reasonably possible. They should be taken 
as a conjecture, and they suggest performing experimental evaluation of the 
exponents for a check of the ideas of the present paper. 



What do we imagine to happen when the equations are changed by enlarging 
the cut off (in the velocity as well as in the displacement variables) or by 
changing the forcing? Suppose that the cut off is already so large that adding 
one extra mode does not really affect the qualitative and quantitative features 
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of the motion. Then adding one mode, i.e. increasing the total dimension of 
the system by 2, should add one pair of viscous exponents at the end of the 
spectrum respectively equal to and — x, as drawn in the third graph in the 
figure. While changing the forcing should, from time to time as the forcing 
changes, change the category of some exponents. Namely the simplest picture 
would be that one of the vanishing slow GED exponents "becomes" positive 
and one of the viscous "becomes" inertial and paired with it; a symmetric 
evolution should take place with the displacement exponents. Or vice-versa. 
The attractor changes dimension by 2 units, sec [BG],[BGG], at each of such 
events. 

In order that the latter picture be possible one needs that at a transition the 
viscous spectrum bottom consists of a pair of a > displacement exponent and 
a negative viscous exponent coinciding with the inertial exponents top pair; 
and at the same should happen for the bottom pair of the inertial and the top 
pair of the slow spectra. 

The case of periodic boundary conditions does not fit in the above analysis 
of the pairing rule because on the torus there is no way of forcing the system 
with a locally conservative but globally non conservative force field with 
average. Nevertheless some kind of pairing might still occurr under simple non 
conservative forcing acting only on some large scale modes, see §5. 

It has been pointed out to me by F. Bonetto that consistency of the above 
picture requires that the sum of the displacement exponents be exactly 0: the 
two of us have indeed been able to verify that this property is formally exactly 
verified in the ED equations. And this led to a correction of the graphs drawn 
in the figure above that I had originally drawn without taking such property 
into account. We shall come back on this point in a future study. Note also 
that the fact that the sum of the displacement exponents vanishes provides 
a natural test that the truncation that one is using is actually large enough 
for having reached cut-off independence of the asymptotic properties ot the 
motions: this happens at the cut off value where the sum of the displacements 
exponents vanishes: further addition of modes only makes longer the flat part 
in the third graph of the figure above. 

§5 The NS and GNS equations. Extension of the pairing rule. 

We turn to the NS and GNS equations, whose interest is far greater than the 
just studied ED or GED equations. 

One is tempted to say that the scenario should be the same. However the 
pairing rule analysis, which is essential for the physical interpretation of the 
results, is no longer naively possible, not even at a heuristic level. 

A pairing rule, first pointed out in special non constant friction cases in [EM] , 
p. 281, has been proved only in the case of systems subject to special gaussian 
constraints, see [DM], but it has apparently a much wider validity, see [EM], 
[ECM1], [SEM], [DPH] and it is likely to hold also in the cases GNS and NS, 
at least in some sense. 

Also the argument in [Dr] implies the existence of pairing in systems that are 
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obtained from hamiltonian systems by adding to them an irreversible constant 
friction term proportional to the momenta in a system of canonical coordinates. 
And the argument in [DM] is restricted to "isokinetic" constraints precisely be- 
cause they are reversible constraints that are obtained by adding to a hamilto- 
nian system a suitable force proportional to the canonical momenta. Since this 
is an essential feature for the validity of (3.2) the latter becomes doubtful in the 
cases (that include NS and GNS equations) in which the friction or thermostat 
forces are proportional to the canonical momenta via a matrix C which is not 
the identity (it is the laplacian in the case of the NS or GNS equations). 3 
In such cases one could envisage that (3.2) is replaced, in the GNS equations 
case, by a relation like: 

^t^L = (f3) Cj (5.1) 

where (/?) is the n v ,gns average (in the case of the NS equations one would write 
v instead of (/?)); and Cj is some suitable function of j, that might be related 
to the spectrum of the matrix C. However attempting at establishing such a 
connection would lead to too many too detailed assumptions at this stage and 
one would like not to rely on them. And from the proofs in [Dr] , [DM] it seems 
unlikely that a pairing rule can hold in a strong form, i.e. that (5.1) holds for 
the local exponents if {(3) is replaced by (3. 

We therefore define Cj by the (5.1) without linking Cj to the matrix C . However 
we shall suppose that (5.1) holds in a "almost local" form in the sense that 
on a rapid time scale (5.1) becomes true also for the local exponents. This 
means that, up to an error that tends to zero very quickly with the time r, the 
logarithms of the eigenvalues of the matrix ( J^(x) J t (x)) 1 / 2t , with J T (x) being 
the jacobian matrix for the evolution operator V T at x, divided by r verify 
i(A+ + Xj) = Cjf3 T (x) with P T (x) denoting the average i J^ 2 f3(V t x) dt. 

This property, together with the Axiom C assumption, suffices to extend, in 
a suitable form, the validity of the predictions (i.e. conjectures) discussed in 
this section for the ED and GED equations to the case of the NS and GNS 
equations as follows. 

We first remark that the really relevant feature of the pairing rule, as far as 
the applications in [BGG] and above are concerned, is not the constancy of 
the pairing but, rather, the fact that some kind of pairing takes place on a fast 
enough time scale, secondly we assume that this is actually the case for the 
GNS equations. These are the main remark and the main assumption on which 
we base the analysis of the fluctuation theorem predictions in the NS and GNS 
equations. 

If a relation (5.1) holds the constants Cj will have to add up to the sum of the 
Lyapunov exponents that can be derived as the average value of a: this means 
that (cr)+ = 2D K {P)r ) .,gns = (X^ = f c J )(f3) n .gns, up to terms ncglegible as the 
Reynolds number tends to oo, see (1.10) and comments preceding it. 



3 I owe to F. Bonetto the clarification of this point. 
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Furthermore let / be the set of P inertial pairs (i.e. of pairs of Lyapunov ex- 
ponents XJ with one positive element) and suppose that the (5.1) becomes 
valid on a sufficiently fast time scale, then the values of (a) + and (uo)+ would 
have ratio (see [BGG]) (J2jei c j)/(J2j c j) so that: 

(<7 ) = = (E^)w^ = (E c > d = Tv ( 5 - 2 ) 

^3 3 jei jei 

having used the conjecture NS of §1 equating (/?),,, gns to v. 
Then if a local time reversal exists on the attractor (i.e. if the geometric 
Axiom C is assumed as well, [BG], for the dynamics generated by the GNS 
equations) the fluctuations of the observable a will have a "free energy" (or 
a "generalized sum of Lyapunov exponents" to adhere to the terminology in 
[FP], [BJPV]) ((p), in the sense of (4.1), with an odd part pPv, with P defined 
in (5.2). This is a property whose validity can be conceivably tested in, real 
or numerical, moderately turbulent systems. At least the linearity in p of 
C(p) — C( — P) should be observable. 

Note also that, in all cases, the pairing rule is trivially valid in the case of no 
forcing: in fact the equivalence criterion in the conjecture in §2 requires that 
in absence of forcing one has to take r\ = or e = 0: i.e. the stationary state 
is, in that case, the trivial (non chaotic) flow u = , 6 = const. 

The assumption that the forcing be locally conservative has not been used and 
disappears togheter with the constancy of the pairing: the above more general 
pairing hypothesis (see (5.1) and the comment following it) is more "flexible" 
and does not require the special hypothesis of local conservativity of the the 
forcing. 

§ 6 Relation between the NS and ED equations. The barometric formula. 
Finally we discuss another main point of our analysis. 

In reference [Ga6] the argument leading to the conjecture NS above can be 
interpreted as saying that NS and ED are also in some sense equivalent. 

The argument is based on the constancy of the dissipation rate e in a stationary 
flow at high Reynolds number and on the microscopic reversibility. In some 
sense the GNS equations emerge as even more natural than the NS equations. 

A criticism can be raised, however. In fact one can argue that the energy is 
also constant in a stationary state and one could develop the argument in [Ga6] 
to imply that the GED equations are also a good model for a fluid motion. 

Since clearly one should not expect NS and ED to be equivalent this looks at 
first as an unsoluble logical contradiction. Which can furthermore be conceiv- 
ably easily checked to occurr. 

However on further thought the contradiction can be resolved and one should 
think that all what has been deduced is that there should be a relation between 
the stationary states of ED (or GED cquivalcntly) and of NS (or GNS). The 
relation to which I think is the kind of relation that one also finds in equilibrium 
statistical mechanics in gases in a strongly varying external field of intensity g, 
like the gravity field. 



1G 



Locally a gas in a field looks just like a homogeneous gas in equilibrium, but 
globally over a length scale H over which the external potential really changes 
(fimgH ~ 1, if (3 is the inverse temperature and m the particles mass) one will 
see that pressure and the density are not constants and one gets the barometric 
formula, see [MP]. 

Likewise we can expect that the stationary states of ED (or equivalently of 
GED) are also "locally" equivalent to stationary states for NS (or GNS): in 
the sense that if we only look at observables depending on field components 
u k with modes k on a certain scale | k | ~ k whose size depends on the dis- 
sipation then we should see essentially no difference. The precise relation that 
determines k will be called barometric formula: it should be easy to determine 
the formula on the basis of dimensional considerations. Locality is here to be 
interpreted in momentum space rather than in coordinate space. 

The determination of the barometric formula amounts essentially at a devel- 
opment of the analogous of the Kolmogorov theory for the ED equations. 

We now attempt at a development of such theory, in the no obstacles case 
for simplicity, on the basis of a few assumptions that we did not subject to as 
much criticism as they certainly deserve: a criticism that will be undertaken in a 
separate publication. We follow closely the ideas (and imitate the assumptions) 
of the exposition of Kolmogorov's theory in [LL] . We set p = 1. 

It seems reasonable to suppose that in the ED case the stationary distribution 
equipartitions the energy among the modes, i.e. (| 7 | 2 ) = -f 2 for all k in the 

"inertial range" L^ 1 <C \k \ -C k x where k x is the "Kolmogorov scale", to be 
determined below. Hence 7 2 (fc x L) 3 = e. 

Then a velocity variation characteristic of the momentum scale k is v 2 — 
k |e[«/2 k] 7 fe ) 2 ) an d 7 assuming statistical independence of the distribution 
of the various 7 , we get v 2 — (kL) 3 j 2 up to a constant factor. 

The scale k x can be determined because it has to be a momentum scale that 
only depends on e and \ so that k x = (p-) 1 / 2 up to a constant factor: this is 
a momentum scale analogous to the Kolmogorov scale k v of the NS equations, 
which is k v = (^) 1/4 . 

For purposes of comparison we note that the quantity called e in the Kol- 
mogorov's theory ("K41-theory"), see [LL], corresponds to of the present 
paper. 

In this case the energy distribution (i.e. the amount K(k)dk of energy per 
unit volume and between k and k + dk) is K(k) = ff p-, for k < k x : very 
different from the Kolmogorov's fc~ 5//3 law. 

In the K-41 theory a key role is plaid by the quantity v^k which is identical 
to r\v for all k <C n -C k v . Therefore we compute the value of v\k in our case 
and we find: 




(6.1) 
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and we see that the quantity v^.k does depend on k in the ED case. Given k the 
SRB statistics for the ED equations driven with a total energy e gives to this 
quantity the same value that it has in the SRB statistics for the NS equation 
driven with a total vorticity n if: 



provided (of course) k is greater than the Kolmogorov scales k Vl k x . 
The "barometric formula" is then the statement of equivalence between NS 
and ED on scale k, i.e. if one only looks at field properties depending on 7 

for i/t < I k\ < k, if (6.2) holds and k 3> k v , k x . 

If we look at a different scale «/ = 2™k for some (large) n then we can expect 
equivalence between ED (or GED) and NS (or GNS) but the pairs e, r\ should 
now be such that (6.2) holds on the new scale: the analogy with the usual 
barometric formula for the Boltzmann Gibbs distribution in the gravity field 
justifies the name given to (6.2). We see that r\v plays the role of the gravity, 
ev plays the role of the chemical potential and n/k x plays the role of the height. 

The above analysis seems to be fully consistent with the numerical results in 
[ZJ] who first proposed, in a different context, a picture very close to the one 
developed here. 

It is clear that this point of view has several consequences: for instance in 
particular it tells us that that the shape of the pairing curve in (5.1) cannot be 
arbitrary {i.e. (P)cj <~ v kj if the modes are ordered in increasing order). This 
is a point on which I hope to return in a later analysis. 

Also: the equivalence between NS and ED on a given momentum scale makes 
more interesting the ideas in [Ga4] and a test of the Onsager reciprocity derived 
in the latter paper seems now quite feasible and seems also to have conseqiences 
for the real NS equations. 

A further remark is that although (6.2) depends on the validity of the K-41 
theory and of the corresponding theory for ED equations the barometric for- 
mula can be developed independently of such theories: hence any modification 
of the K-41 theory (and of the corresponding theory for ED) will lead to a 
barometric formula, with a relation between e, 77, k, possibly more complicated 
than (6.2). 
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Appendix Al: The hamiltonian formalism for Euler equations and Dressier 's 
theorem for ED. 




(6.2) 
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To check the applicability of the results on pairing of [Dr] to the ED equations 
we must check that the equations can be written, in canonical coordinates for 
some hamiltonian function H, in the form: 

£ =d p _H 

P=- 9 g H+ F-xp 

where F is such that d qj Fi = diFj without being necessarily F_ = — dV 
for some globally defined V (the latter would be a trivial case). The labels 
for the components of q are x,i with x £ f2 and i — 1,2,3. The partial 
derivatives are, correspondingly, functional derivatives; we shall ignore this 
because a "formally proper" analysis is easy and leads to the same results. By 
" formal" we do not mean rigorous, but only rigorous if the functions we consider 
have suitably strong smoothness properties: a fully rigorous treatment is of 
course impossible for want of reasonable existence, uniqueness and regularity 
theorems for the Euler equations or the Navier Stkes equations in 3 dimensions. 
Consider first the Euler equations. They can be derived from the lagrangian: 

Co(8,S) = £ [ 5 2 dx (A1.2) 

(p = density) defined on the space V of the diffeomorphisms x — > S(x) of the 
box f2, by imposing the ideal constraint: 

det J= det-^(x) = dSx A d& 2 • dS 3 = 1 (Al.S) 

ox 

In fact, if Q is a Lagrange multiplier, the stationarity condition for: 

£{L I) = 7; J l 2 dx+ I (}(r)dv{ J(S )(r)<U ( 11.4) 
leads to: 

pS = -(det J) J^dQ = -det J -^=- -d Q {Al.S) 

so that setting u(S(x)) = S(x),p(S) = Q{x) if 5 = S(x), we see that: 

^L = --d P (A1.6) 

dt p- F x ' 

which are the Euler equations. And the multiplier Q(x ) can be computed as: 
Q(x)=p(5) = -[A- 1 (du -d u)] 5 _ (A1.7) 

where the functions in square brackets are regarded as functions of the variable 
8 and the differential operators also differentiate over such variable. After the 
computation the variable 8 has to be set equal to S(x). 



9maggiol996 



19 



Therefore by using the lagrangian: 



Ci{Li) = J {^f^-[A-\du -8 «)] 4(£) (detJ(5)| £ -l))ds (A1.8) 

we generate lagrangian equations for which the " surface" E of the incompress- 
ible diffeomorphisms in the space D is invariant, these are the diffeomorphisms 
x — > <i(x) such that J(<5) = d5\ A 862- 863 = 1 at every point 2; G 0. 
Then E is invariant in the sense that the solution to the lagrangian equations 
with initial data "on E", i.e. such that 5 G E and 9 • j (x) = 0, evolve 
remaining " on E" . 

The hamiltonian for the lagrangian (A1.8) is obtained by computing the canon- 
ical momentum p(x) and the hamiltonian as: 

- ( -> (AL9) 
H(p,q)=±(G(q)p,p) 

where G( q ) is a suitable quadratic form that can be read directly from (A1.8) 
(but it has a somewhat involved expression of no interest for us), and the 
. . . (that can also be read from (Al.8)) are terms that vanish if 5 G E and 
d ■ 5 — 0, i.e. they vanish on the incompressible motions. 
The above is well known and shows that the Euler flow can be interpreted as 
a gedodesic flow on the surface E of the incompressible diffeomorphisms of the 
box enclosing the fluid, see appendix 2 in [A]. 

Modifying the Euler equations by the addition of a force / ( x ) such that 
locally f (x ) = — d&(x ) means modifying the equations into: 

^ = -d P -d^ (ALIO) 
which can be derived from a lagrangian: 

jCf(S,S))=jCi(6,6))- J $(5(x))dx (ALU) 
which leads to the equations: 

u(5(x)) = -- p d 5 _p(S(x)) +d & *(8(x)) (ALU) 



Hence the ED equations have the form: 



q =dpH 

p =-d q H - F-XP 



(A1.13) 



at least as far as the motions which have an incompcssiblc initial datum are 
concerned. This is true because the ED equations which have an incompressible 
initial datum evolve it by keeping it incompressible. 
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The Lyapunov exponents of the equation (A1.13) verify the pairing rule by 
the analysis in [Dr] . However the pairing takes place in the full phase space of 
the diffeomorphisms of ft, including the incompressble ones. 

It is not difficult to see, by using that the constraint to stay on the surface £ is 
holonomous, that one can find canonical coordinates 7r , n , 7r , k - 1 describing 
the motions on S or, respectively, transversally to it. And the equations for 
7r k - 1 arc, near X and for n - 1 small, ir_ = — XTL 1 ' an d ^ — TL^ so that 
the corresponding Lyapunov exponents are trivially paired in pairs 0, — \ with 
pairing sum — ^. 

Since we have seen above that all the exponents are paired at the level ^ this 
means that all the physically intereting exponents (relative to the incompressible 
motions, i.e. relative to the w , k coordinates) are also paired at the same level, 
as claimed in §4. 
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